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Abstract 

We propose a new approach for describing the irreversible behavior observed in small quantum 
systems. It is based on the fact that equilibrium thermodynamics may reveal the existence of 
an underlying dynamics. This is true in the algebraic approach of quantum mechanics via the 
Tomita-Takesaki theorem. A similar result is obtained if we start from the path integral expression 
of the partition function, then an equation of motion is introduced. It corresponds to a continuous 
diffusion process but retaining some aspects of the uncertainty relations. We discuss the relation 
between this equation and the Schrodinger equation. Although the equation of motion is time- 
irreversible it can be used to describe exactly all the equilibrium properties of a thermodynamic 
system. The equilibrium is defined in such a way that there is an exact compensation between the 
energy spent, in average, on the paths and the energy put in the system from outside. From this 
equilibrium condition we may give a meaning to the thermal time in path integral formalism. In 
a second step we use the equation of motion to describe irreversible dynamics. We show that the 
equation of motion is equivalent to a quantum Smoluchovski equation. The relation between our 
Smoluchovski equation and similar ones used in the literature is emphasized. From the comparison 
between the relaxation time that the system needs to reach its equilibrium state and the thermal 
time we show in which condition the concept of thermostat can be relevant. A standard bistable 
model is investigated. The chemical rate is calculated as a function of time, it appears to be a 
non-monotonic function of time. In very particular conditions the stationary value of chemical rate 
can be identified with the Kramers result. From all these results it appears that our equation of 
motion or its Smoluchovski version may be considered as a realistic starting point for describing 
quantum dissipation in small systems. 

PACS number:03.65Ca, 05.30-d, 05.70-a, 47.53+n . 
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I. INTRODUCTION 



To describe the irreversible evolution of a system the first fundamental element we have 
at our disposal is the second law of thermodynamics from which the so-called arrow of 
time is introduced. This law asserts the existence of a state function, the entropy, that is a 
non- decreasing function of time for any closed system (-). To derive the second law from sta- 
tistical mechanics we have to explain how a macroscopic system may exhibit an irreversible 
behavior when the microscopic description is based on time reversible equations (see for 
instance (-)). Boltzmann tried to solve this problem via its derivation of the H — theorem 
that generated an extensive literature (-). To explain the irreversible behavior of chemical 
reactions but staying in classical physics Kramers (-) suggested to use the Fokker-Planck 
equation that results from general aspects of the theory of random processes In this 
equation position and momentum of particles appear on the same footing. In some cases 
it is possible to eliminate the momentum and to get a Smoluchovski equation. Different 
regime can be observed depending on the value of the friction coefficient (-). 
In modern approaches of chemical kinetics the starting point is a quantum description. Then 
we have to deal with the following question: how to describe quantum dissipation if we use 
the Schrodinger equation that is time reversible in the Wigner sense. Quantum mechan- 
ics is based on a hamiltonian formalism which can not describe dissipation processes (-). 
The system-|-reservoir methods give a route for introducing irreversibility (see for instance 
(-)) and they have permitted a large improvement in our understanding of the dynamics in 
dense medium. This kind of approach has been extensively used to describe, for instance, 
the tunneling of electrons through a dissipative barrier a very important problem in simple 
electronic devices (-). 

In the system+reservoir approaches the system under investigation is considered as the 
small open part of a large system. The dissipation arises from the energy transfer from 
the small system to its large environment. In such approaches the dynamics is described 
via a quantum mechanical Langevin equation (see for instance (-)) for the relevant opera- 
tors of the reduced system or via a generalized quantum master equation for the reduced 
density matrix. Frequently the master equation is a Fokker-Planck type equation (-). The 
system-|-reservoir approach gives a correct description provided the system relaxation times 
are larger than t = (ih where /3 = l/Zc^T and h is the Planck constant (-). In absence of 
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reservoir the system is described by the Schrodinger equation. 

At this point it is important to underhne the difference between a reservoir considered in the 
system+reservoir methods and the usual meaning of resevoir in standard statistical mechan- 
ics. In the last case the reservoir is an infinite system without any specific property, its main 
role is to fix some external measurable parameters like temperature or chemical potential. 
If the reservoir is used to fix the temperature T it is called thermostat. In system+reservoir 
approaches specific interactions exist between the system and its reservoir and by fixing the 
properties of these interactions we may reproduce some well known properties of the small 
system, for instance the existence of a brownian motion in the classical limit (-). In order 
to eliminate the reservoir an average over the reservoir variables is performed by using the 
canonical density matrix assuming that the reservoir is in contact with a thermostat. Here- 
after we associate to "reservoir" the meaning it has in the sytem+reservoir methods while 
"thermostat" corresponds to the standard concept in statistical mechanics. 
In what follows we want to investigate the irreversible behavior of small systems isolated 
from any reservoir but eventually in contact with a thermostat. In standard chemical kinetics 
we have in mind the description of a reaction such as 

A + BC ^AB + C (1) 

involving three chemical species A, B, C. Our main goal is to find an equation for describing 
the irreversible behavior of such systems. Using the concept of reaction coordinate a 
traditional analysis of ([1]) leads to consider the properties of a system in presence of a 
bistable potential (-). We will retain this kind of approach here and we will analyze the 
case of a particle inserted in a box in which there is a bistable potential. Such a system has 
been investigated in the seminal paper of Kramers ('^) in order to mimic a nuclear fission 
process. More generally molecular double minimum processes have been the subject of 
many experimental and theoretical investigations, in the attempt of a better understanding 
of elementary chemical reactions. 

This paper is organized as follows. In Section 2 we focus on systems at the thermodynamic 
equilibrium and we show that for such systems it exists an underlying dynamics. This idea 
has been first developped in the so-called algebraic approach of quantum physics in relation 
with the Tomita -Takesaki theorem (-), here a similar idea will be extended to the path 
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integral formalism following an idea suggested by Feynman (-); an equation of motion is 
derived. We analyze the physical content of this equation; in particular we discuss the 
relation between this equation and the one obtained by introducing an imaginary time in the 
Schrodinger equation. The equation of motion we consider is associated with a continuous 
diffusion process but retaining some aspects of the uncertainty relations. From the entropy 
a thermodynamic equilibrium condition is defined and we give a meaning to the thermal 
time in path integral formalism. Our main assumption is that this equation of motion 
is the general equation governing the evolution of a particle in an external potential. In 
Section 3 we show that the equation of motion is equivalent to a quantum Smoluchovski 
equation; its solution is a density of probability, in which a Boltzmann distribution is mixed 
with a quantity associated with the path integral formalism. We may give a meaning to 
this probability and we discuss the possibility of coupling our system with a thermostat. 
In Section 4 we use the Smoluchovski equation to describe the irreversible evolution of a 
particle injected in a box in which there is a bistable potential. This standard model has 
been investigated from a general point of view in (-), (— ) and applied, for instance, to 
describe nuclear fission processes in (— ). For this model we investigate the chemical rate 
and its time dependence. Conclusions are given in Section 5. 

II. DYNAMICS IN THE THERMODYNAMIC EQUILIBRIUM 

Frequently, when a new field of investigation is developing different equivalent approaches 
are proposed. By focusing on different starting points they may reveal complementary 
aspects of the theory. This can be illustrated in the the case of the density matrix. As 
summarized below from the recent developpments of the matrix mechanics elaborated by 
Heisenberg it has been concluded that a given dynamics is associated with the existence 
of the density matrix. We show that an underlying dynamics can be also associated with 
the path integral formalism, it is described by an equation of motion. A large part of this 
Section will be devoted to an analysis of the physical content of this equation. 
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A. Thermodynamic equilibrium in algebraic quantum mechanics 



The matrix mechanics has given rise to the so-called algebraic approach of quantum me- 
chanics (— ) in which the main point is the existence of a C* algebra of non-commutative 
operators. An important result of this approach is the Tomita-Takesaki theorem (for a 
recent review see (^^)) that establishes the existence of a one parameter group of automor- 
phisms that leaves the algebra globally invariant. A consequence of the Tomita-Takesaki 
theorem is the existence of a relation 

{B,{at+,A)) = {{atA),B) (2) 

in which (, ) means the inner product and at is the so-called modular group on the algebra. 
At this level it is impossible to claim that the parameter t in ([2]) is a time. However we 
may compare ([2]) with the KMS condition (— ) established by Kubo and Martin and 
Schwinger (— ). Haag and his coworkers i}'^) postulated that the KMS condition is the 
correct definition of thermal equilibrium for infinite dimensional quantum systems. The 
KMS condition shows that the correlation function between two variables A and B noted 
< {jtA)B > is analytic in the strip < Im{t) < (5% and we have 

< {^tA)B >=< Bi^t+i^nA) > (3) 

where jt introduces the time translation group defined according to 

itH ^ itH 
'jtA = exp — ^ exp — — (4) 

in which H is the hamiltonian operator. From the comparison of (|2]) and ([3]) it was con- 
cluded that the Tomita-Takesaki theorem is equivalent to the time evolution of the bounded 
operators generated by the Hamilton provided the time is measured in units Ph (^^). 
These results show that the thermal equilibrium characterized by an algebra of operators 
contains an underlying dynamics; these results depend on the form of the density matrix 
but they do not refer explicitly to the existence of the Schrodinger equation. 



B. Thermodynamic equilibrium in the path integral formalism 



Starting from the canonical density matrix 



exp — (3H 
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where a normalization constant. Feynman (-) derived a path integral formalism for the 
partition function Z . In order to do that it was needed i) to start from the Schrodinger 
equation from which we must calculate stationary states and energy eigenvalues, ii) to de- 
velop some arguments justifying the use of the canonical form of the density matrix and iii) 
to introduce some mathematical tricks. Then the partition function can be written as 



in which Dx{t) is the measure for the functional integral and y4o[x(t); 0, r] represents the 
quantity 



calculated on loops i.e. on trajectories on which we have x{t = 0) = x{t = t) = x; for each 
loop the dummy variable, t, varies from to r; t and r have the dimension of a time. In 
([7]) it is assumed that the external potential u{x{t)) does not depend explicitly on t. The 
integral (|6]) is calculated on the all values of x taken in the volume of the system under 
consideration. The expression of Z given in is just a mathematical transformation of 
the standard expression of Z = Tr(exp[— provided we use r = ph. Of course in such 
a formal derivation we cannot claim that t is a time having a physical meaning i.e. that 
Ao[x{u); 0, r] is an action associated with loops in ordinary space time. 
After deriving ([6]) Feynman (-) wrote a paragraph entitled " Remarks on methods of deriva- 
tion" on which he suggested that it should be possible to derive the expression of Z directly 
from the description of the motion how he did for the quantum mechanical amplitude in 
quantum mechanics. Why to search such a short cut that should avoid the wave function 
and the energy levels ? The Feynman's answer is the following "in doing that a deeper 
understanding of physical processes might result or possibly more powerful methods of sta- 
tistical mechanics might be evolved". Hereafter, to be short we will mention these remarks 
as the Feynman's conjecture {FC). From the FC it is suggested that Z = Tr{exp[—PH]) 
contains the actual dynamics of the system provided we use its path integral version. This 
is in the same spirit of what has been done in the algebraic approach of quantum mechanics. 
In order to find the dynamics associate with Z we define the quantity 




(6) 









q{to,Xo]t,x) = J Dx{t) exp --A[x{t);to,t] 



(8) 
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in which A[x{t);to,t] is defined in ([7]) but not restricted to close paths. From q(to,Xo;t,x) 
and a well behaved function (f)o{x) defined for [t = to) we may form a real- valued function 
(j)(t, x) according to 

x) = J 0o(z/)g(^o, y, t, x)dy (9) 
By using the Feynman-Kac formula, we can see that x) is the solution of the equation 

h 1 

- d(p{t, x) /dt + — A^0(t, x) - -u{x{t))(f){t, x) = (10) 

that verifies the initial condition 0(0, x) = (poiv) i-^- (l{to,xo',t,x) is the fundamental 
solution of ffTOl) in which means the laplacian operator taken at the point x. 
Hereafter we consider the equation of motion (ITOl) as the starting point of our approach. 
We cannot proove that ffTOj) is the actual equation of motion but we can try to justify it 
i) by analyzing consequences that we can derive from it and ii) by showing that it is in 
agreement with general trends in modern physics. 

C. Physical meaning of the equation of motion 

From the canonical form of the density matrix the previous approaches lead to introduce 
the usual evolution operator but with a rescaled time in the case of the algebraic approach or 
to the equation (fTOj) in the case of the path integral formalism. In principle these two results 
do not refer explicitly to the Schrodinger equation. Indeed, in the spirit of the FC it must 
be so and ( ITOl) must be justified by itself from the underlying dynamics. Nevertheless we 
observe that (fTOl) corresponds formally to the Schrodinger equation provided we introduce 
an imaginary-time, this suggests that a link must exist between ( ITOl) and the Schrodinger 
equation. 

1. Relation with the Schrodinger equation 

The imaginary-time approach is based on the observation (^) that the partition function 
is formally identical to the quantum mechanical amplitude provided the time t is replaced 
by it. By this purely formal procedure at least one important question remains unsolved: 
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why the solution of (fTOj) that is a real valued function has a physical meaning while only 
the product of the wave function by its complex conjugate has a meaning ? We must search 
a more fundamental relation between (IIUI) and the Schrodinger equation. This is what we 
have done in by focusing on the reversible/irreversible behavior of these two equations. 
To have a self contained paper, hereafter we shortly recall the physical content of this work; 
the main point being how to complete the dynamics given by (fTOj) in order to describe a 
reversible behavior. 

The Schrodinger equation is reversible in the Wigner sense, this means the following. If 
ip{ti, x)) is a wave function, its evolution is given by a unitary operator U{ti, to) relating the 
wave function taken at time ti and the one at time to according to ip{ti, x) = U {ti,tQ)ilj{to, x). 
If at time ti we consider 'ip*{ti, x) the complex conjugate of ip{ti, x)) then by using the same 
evolution operator and its unitary character we can show that for a time interval (ti — to) 
counted after ti we get a function 4'*(to, x) which is the complex conjugate of ip(to, x). Thus, 
to get a Schrodinger equation from ( fTOl) we must introduce a second equation describing a 
backward motion but with the same fundamental solution as ( ITOl) . its solution will be noted 
(j){t, x). A fundamental result due to Kolmogorov (see (— ) and the references quoted therein) 
is to show that the equations describing forward and backward motions are in duality and 
from this it is possible to define a Markov process for which x)(j){t, x) is a density of prob- 
ability. The two real valued functions (j){t, x) and 0(t, x) can be combined into a complex 
valued function ilj{t,x) verifying a Schrodinger equation (■^) and the Born interpretation 
of the wavefunction becomes obvious. Thus we see how from physical arguments we can 
recover a Schrodinger description from ( ITOl) . 

Now we have to deal with the following problems: what to say about the dynamics repre- 
sented by f|TO|) ? can we justify this equation ? can we associate to it an acceptable physics 
? These problems will be considered in the next subsection but detailed arguments already 
developped in (— ), (— ) will be not reported here. 

2. The dynamics associated with the equation of motion 

There are many indications showing that spacetime may be discrete rather continuous 
(see for instance (^^) and the references quoted therein). The choice between a discrete and 
a continuous version of the spacetime structure has already been analysed by Riemann in 
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the classical world and more than 30 years ago Feynman presented some doubts concerning 
the continuum nature of spacetime in the quantum domain (^^). Today it is well accepted 
that the conventional notions of space and time break down at the Planck scale where new 
uncertainty relations have to be introduced. A discrete spacetime means that any length is 
built up from a finite number of the elementary length, 5x, and any time interval results 
from a series of individual ticks of duration 5t. To give a structure to the spacetime we 
have to introduce a relation between 5x and 5t this relation determines the scale at which 
we want to describe the physical world. At the Planck scale the new uncertainty relations 
play this role and it is stablished that there is a minimum for 5x and 5t depending on the 
velocity of light and on the gravitation constant. In the pre-relativistic domain in which we 
are concerned the only one universal constant that we have is h and for a mass m the only 
one relation that we can introduce between 5x and 5t is = h/m. This relation is a form 
of the Heisenberg uncertainty relations since we have immediatly 6x6p = h and 6t.6E = h/2 
provided we use 6p = ■m{5x/5t) and 6E = (l/2)m(^)^. 

Although this is probably not needed we assume that the spacetime points are located on 
the sites of a regular lattice, as in the chessboard problem investigated in (■-), the lattice 
constants being 6x and 6t. In absence of external field we assume that the motion is as 
simple as possible. A particle may jump, at random, from one site to one of its nearest 
neighbours. This corresponds to a random walk. In absence of extra conditions fixing the 
values of 6x and 6t and we may assume that 6x and 6t tend to zero however in this limit 
the ratio {6xy/{St) must remain finite and his value is h/m. The equation of motion is 
equivalent to a continuous diffusion process in an external field and the diffusion coefficient 
D = h/2m is related to the existence of uncertainty relations. 

3. Properties of the equation of motion 

The solutions of ffTOl) are real valued functions <^(t,a;). However since ffTOl) is not a 
Chapman- Kolmogorov type equation due to the presence of the external potential the inte- 
gration of (j){t, x) over a finite volume is not a time independent quantity and, as a conse- 
quence, (f){t, x) can not be normalized. The function (j){t, x) is a weighted sum of all the paths 
arriving at the point x and at time t when the initial condition (j)o{y) is given. The equa- 
tion (fTOl) is time- irreversible and from it the associated dynamics corresponds to a positive 
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semi-group for which we can show (^°) that q{to,xo', t,x) verifies the law of composition 



q{ti, xu t2, X2) = / dx3q{ti,xu ts, X3)g(t3, X3; ^2, X2) 



(11) 



provided ti < ta < ^2 and therefore it can be used to describe the transitions in space-time. 
In (— ) we have shown that all the averaged quantities calculated with ffTOj) on closed loops 
corresponding to an equilibrium situation have a clear physical meaning. 

4- Dynamic point of view of the thermodynamic equilibrium 

We have seen that a simple dynamics is associated with ffTOj) . However, in the spirit of 
the FC, we have to explain how from this time irreversible equation we may describe the 
thermodynamic equilibrium. This is the main goal of this subsection in which we will define 
the meaning of the dynamic equilibrium and the physical sense of the thermal time in the 
path-integral formalism. In order to do that we will start from a definition of the entropy 
since we know that the entropy is the corner stone from which we may describe the thermal 
equilibrium (— ) . 
Let consider the quantity 



in which r is, for the moment, an undetermined free parameter having the dimension of 
a time and U is the external energy put into the system during its preparation, from a 
thermodynamic point of view U corresponds to the internal energy of the system. In (fT2|) the 
quantity tU has the meaning of an external action. If it exists only one possible trajectory 
for which the euclidean action Ao[x{t); 0, r] exactly compensates the external action tU we 
may say that there is no disorder in the system. For real situations it exists a lot of paths for 
which the order of magnitude of Ao[x(t); 0, r] — rU is approximately h, all these trajectories 
contribute to S. Larger is this number of trajectories smaller is the order in spacetime and 
larger is S. The calculation of S requires to start from a point x = x{t = 0) in space and 
to explore during a time interval r all the loops around this point and, finally, to perform 
the same procedure for each value of x in the volume of the sample. It is clear that S 
characterizes the order or the disorder in spacetime, but at this level S is not the thermal 
entropy. 




(12) 
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The quantity 5* defined above depends on two external parameters r and U. From the 
definition f lT2|) we may calculate the derivative In (^^), it has been shown that we have 



in which < uk{x) >path is the regular part of the mean value of the kinetic energy calculated 
over the paths localized around the initial point x and < up{x) >path is a similar quantity 
but associated to the potential energy coming from the external potential. From ( IT3ll we 
may introduce an equilibrium condition 



This condition means that the mean value of the energy calculated on the paths during a 
time interval r is equal to the internal energy needed to create the system, i.e. in average, we 
can not spend more energy than the energy put initially in the system. Such an approach led 
to the concept of thermodynamic time in general relativity ). At the thermal equilibrium 
we know that ^ is a measurable parameter of the system that corresponds to the reverse 
of the temperature and from (fT3l) and (fT4l) we get r = ph. Thus the equilibrium condition 
( |T4l) lead to introduce a particular time scale whic corresponds to the thermal time. Note 
that between t = and t = r no entropy is created in the system as shown in ). When 
this value of r is introduced into the expression of the entropy we recover exactly all the 
thermodynamic results in the path integral formalism (^). 

Clearly r is a characteristic of the equilibrium state, it is not surprising that this unit of 
time also gives the scale of time in the Tomita-Takesaki theorem that it is concerned with 
equilibrium situations. Thus, although the motion on the paths are described by a time 
irreversible equation we can describe a thermal equilibrium via f |T4l) . Note that r does not 
represent the relaxation time that a system taken in a non equilibrium state requires to 
reach its equilibrium state. Such a relaxation time will be investigated in Section4. Another 
example of relaxation towards an equilibrium state has been given in papers devoted to the 
derivation of a if — theorem (—),(")■ 

Another meaning of r can be given from the time-energy uncertainty relation (— ). If t < r 
the quantum fluctuations are larger than fcsT, i.e. the typical value of the thermal energy 
and we can not say that there a well deflned thermodynamics for such short times. This 
gives a quantitative meaning to a very well known idea (-). 



h dS 
ks dU 



r+[U - f dx[< Uk{x) >path + < Up{x) >path] 



dU 



(13) 



U = fdx[< Uk{x) >path + < Up{x) >path] 



(14) 
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5. Our main assumption 



In this Section we have seen that (fTOj) can be considered as describing a continuous 
diffusion process retaining some aspects of the uncertainty relations. Although this equation 
of motion is time-irreversible it allows us to reobtain all the equilibrium properties via the 
equilibrium condition (fT4|) . In previous papers we have seen that (11 01) can be completed by 
a second equation if we want to describe a reversible process, this establishes a link between 
(fTOj) and the Schrodinger equation based on physical arguments. In (— ) we have shown that 
all the quantities calculated at equilibrium on the paths have a physical meaning. Thus f lTOj) 
is in agreement with everything we know concerning the thermal equilibrium, hereafter we 
will assume that flTU]) can be also used to describe non-equilibrium state. First, we will show 
that (|T0|) is also equivalent to a quantum Smoluchoski equation. 



III. A QUANTUM SMOLUCHOVSKI EQUATION 

Due to the presence of the potential u{x) in flTUl) . the solution of this equation is not a 
density of probability. However we may use a transformation introduced by van Kampen 
(— ) to establish a relation between ( ITO!) and a Smoluchovski equation and then to define a 
density of probability. 



A. Derivation of the Smoluchovski equation 

Let consider the quantity 

-6V(x) 

P(t, x) = </)(t, x) exp (15) 

in which a potential V{x) is introduced and {0)~^ is a scale for this potential. 
Since V{x) is assumed to be independent of t we have 

dP{t,x) ^ d^{t,x) ^ -evjx) 

dt dt "^^^ 2 ^ ^ 

By replacing ^''^g*'^'' by its value obtained from ( ITOj) and using simple mathematical trans- 
formations it is possible to get the following Smoluchovski equation 

^^1^ = 2^^^'^^'' " ^V.((V.^r (x))P(t, X)) (17) 
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in which V{x) is the solution of 

1 1 2m 

-AJV{x) - -{V.eV{x)r + j^,U{x) = (18) 

Finally we introduce the quantity Z(x,t) = exp ^^y^, it verifies the following equation 

^A,Z{t, x) - l-uix{t))Z{t, x) = (19) 

which is nothing else than ( ITOj) in the stationary regime. Thus from 0(t, x) we can create 
a probability P{t, x) given by f|T5l) which depends on the number of paths arriving at x via 
x) and on the value of the potential V{x) through exp — ^^y^. This is a non-traditional 
expression for a density of probability and we have to give a meaning to P{t, x). 



B. The meaning of P{t, x) 

To solve an equation like (fTOj) we separate the variables x and t searching a set of solutions 
of the form x) = fn{t)ipn{x). The time-dependent functions are given by /„(t) = 
exp( — •^""•^" t) where Eq is the energy of the fundamental state and En is an eigenvalue of 
the equation 

^A.^x) + {{En - Eo) - u{x)){Mx) = (20) 

which is identical to a stationary Schrodinger equation for a particle in presence of the initial 
external potential u{x). Note that the functions are monotonic decreasing functions 
of time in contrast with the solutions of the Schrodinger equation that should be oscillary 
functions of time. Using the closure relation between the eigenfunctions we can write the 
fundamental solution of (fT7|) as 

P(0, y; t, x) = ^,{x)n^]cpn{x) e^-^^^^t) (21) 

where the sum runs on the all values of n from to infinity and P(0, y\ 0, x) = 6{x — y) as 
a consequence of the closure relation. If f{y) is the initial distribution the we have 

P{t,x) = JP{0,y;t,x)f{y)dy (22) 

If t tends to infinity the limit of P{t,x) is Peq{x) = ^Po{xy = cexp —9V{x). In agreement 
with (-), (— ) and (— ) we can interpret Peq{x) as the Boltzmann equilibrium density of 
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probability to be at the point x where there is an effective external potential V{x), c is a 
normalization constant. More generally we can rewrite P{t,x) as 



Thus, before reaching an equilibrium situation, P{t, x) is the product of a Boltzmann distri- 
bution by a quantity '~f{t,x) associated with the paths. In 7(t,x), for each state n, we first 
compare the number of paths arriving at x via t) to the same number at equilibrium 
given by (j)o{0,x) = fo{x), then we perform a sum over the states weighted by 



If at a given time t there is no path arriving at the point x then P{t, x) vanishes and there 
is no chance to find a particle at this time in this point. From all the arguments developed 
above we may interpret P{t, x) as the density of probability to be at the point x at the time 
t taking into account that the initial distribution is given by f{y)- 

If a Smoluchovski or a Fokker-Planck equations have been derived in the system+reservoir 
methods, these equations are based on ingredients including explicitly some properties of the 
reservoir such as the frequency spectrum for the reservoir oscillators or the friction coefficient 
induced by the reservoir particles. The philosophy in the system+reservoirs methods is to 
fit the system reservoir interactions in order to reproduce, for instance, the brownian motion 
for the system in the classical limit (-). In the present work we have no such a reservoir 
and accordingly no such parameters are present in our Smoluchovski equation. In particular 
there is no temperature in our equation in contrast with what happens in equation based on 
the system+reservoir approach or in the classical approach of Kramers (^). The possibility 
of introducing the temperature will be analyzed in the next Section. 

C. Contact with a thermostat 

The equations f|T7j) and f|T8|) show that only the quantity 6V{x) solution of f|T9l) is relevant 
in the calculation of P{x,t). In (-), (— ) the system is assumed to be in contact with a 
thermostat and 6 = P has been chosen. To analyse the relevance of the thermostat concept 
in our approach we have to compare the system relaxation time and the thermal time r = ph. 
From ( l2Ti) and provided that En is an increasing function of n the relaxation time with which 




(23) 



My) 



)f{y)dy 



(24) 
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P{t,x) will reach its equilibrium value Peq{x) is given by (^Ei-Eq) • 

If ^^^^^^^ << (3h or fc^T << (E'l — Eq) the relaxation of P{t, x) takes place in a regime for 
which the quantum fluctuations are much larger than the thermal ones; in such conditions 
the concept of thermostat is irrelevant. Then we have to consider the system as isolated or 
possibly submitted to a distribution of initial conditions. 

In the reverse situations i.e if ^^^^^^^ >> f3h we may assume that the system in contact 
with a thermostat and we have to take a thermal average of (pTI) . This can be realized by 
introducing in (12T]) for each value of n an extra thermal factor given by 

exp(-/3E„) exp(-/3(E„ - Eo)) exp(-/3(K - Eq)) 



Sexp(-/3(E„) Sexp(-/3(E„ - Eo)) 
in which we have introduced the partition function 



(25) 



Zo = Eexp(-/3(K-Eo)) (26) 
Due to the thermal effects P(0, y; t, x) has to be replaced by [P{0, y; t, x)]th that we can write 



as 



[P(0, y- 1, x)U = ^Mx)n^]M^) exp{-^^^-^{t + ^h)) (27) 

We see from f l27l) that the thermal effects change P{0,y;t,x) in two ways; first P{0,y;t,x) 
is rescaled by the partition function Zq and second the time dependence is shifted by the 
thermal time. 

The main result of this Section is to show that the equation of motion f lTOj) is equivalent to 
a quantum Smoluchovski equation ( IT71) in an external potential given by ( ITSll . The solution 
of this equation, P{t, x), gives the probability to be at x for the time t. In the next Section 
we use this equation to investigate an example. 



IV. DYNAMICS IN A BISTABLE POTENTIAL 



In order to illustrate our approach we study the dynamics of a particle injected in a box 
in which there is a bistable potential. This is a standard model investigated in general in 
(-), (— ) and applied, for instance, to describe the nuclear fission (— ). 

Let be a one dimensional box located in the interval —b<x<b,m the region —a < x < a 
with a < b there is a repulsive barrier of height Ui while in the remaining intervals [—b, —a] 
and [a, b] it exists an attractive potential of magnitude Uq. In addition, located at a; = ±6 
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we put an infinite repulsive barrier. Thus the potential u{x) introduced in ( ITOj) is defined 
by four parameters {Uq, Ui, a, h). Our first task is to calculate the quantity 6V{x). 



A. Properties of 9V{x) 

The functions ^Pn{x) are obtained by the usual methods of quantum mechanics and 9V{x) 
can be calculated from ^q{x). This function is given by ^q{x) = Acoshii'ox if —a < x < a 
and <fo{x) = B sin kQ{h — x))iia>x<h the values of the coefficient have been given in (— ). 
It is easy to see that 9V{x) has a local maximum at a; = given by 9V{Q) = (27ra^ + 2 In a) 
with = (3^) and c = (6 — a). 9V{x) exhibits two minima located at ±xo = ±(a + b)/2 
corresponding to 9V{xo) = —2 In 2. The shape of 9V{x) is reminiscent of the one of u{x) 
but the ratio = (— f;^^— ) is not a simple function of a but it also depends on the spatial 
distribution of u{x) through the parameters a and c. Near x = and x = xq, 9V{x) can be 
expanded according to 

9V{x) = 9V{0) - Klx^ = 9V{0) - 9^{2nnofx^ (28) 

9V{x) = 9V{xo) - kl{x - x^f = 9V{x^) - 9]^{2txuj^)\x - x^f (29) 

showing the oscillatory character of the potential near the extrema; the relation between Kq 
and /co and the spatial frequency VLq and ojq are the same as the ones used by Kramers (-). 



B. The chemical rate 

We can calculate the chemical rate from the definition (— ) 
in which 

P{t) = f Pe,(x)Sc„(^^^)rfx (31) 

the integration is performed on the domain —b<x<0. To have a simple model we assume 
that the initial distribution is centered at x = —{b + a)/2 and its extension, a, is very small 
in comparison with (6 — a). 

In the stationary regime i.e. for very large values of t the time-dependent part of P{t) is 
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restricted to -Pi(t) = / Peq{x)ci{^^^^)dx and the rate constant is given by ksta = 
which is nothing else than the reverse of the relaxation time introduced in SectionSC . By 
using the values of {Ei — Eq) given in (— ) and the value of 6*^(0) we obtain 

ksta = -^^o-fco exp -eV{0) = 2TT{—)eno.uo exp -9V{id) (32) 

This result looks like the one obtained by Kramers since ksta is determined by the prod- 
uct VLq.ojq of the frequecies of the potential near extrema (see and (12^ ) multiplied by 
exp— ^^V(0) that we can considered as the probability to cross the barrier. Indeed ( !32|) 
becomes identical to the Kramers result (-) if we take 6 = P and and choose the friction co- 
efficient 7] to be such as ^ = -^9 of course in this case the Smoluchovski equation used here 
becomes identical to the one investigated by Kramers (equation (11) in (-)). Nevertheless, 
as already noted in (^^), this result is not trivial because Kramers used a phenomenological 
trick (-), (-) instead of solving the Smoluchovski equation. 

However it is important to note that fl32|) is not identical to the one of Kramers in gen- 
eral. In the equation of Kramers as well as in the Fokker-Planck equations deduce from 
the system+reservoir approaches, the temperature appears at the level of the Smoluchovski 
equation. In our approach ( l32l) gives the rate constant for an isolated system. The couplig 
with a thermostat is relevant provided we have ^^^^^^^^ » f3h as discussed in Section 3C 
then we may replace ksta by its thermal average given by 

_^ exp{-P{E,-Eo)) ^ Er-Eo 

[Kstalth - [ J ^ [<i<i) 

Our result is richer than the one of Kramers since we have an exact solution of the Smolu- 
chovski equation and consequently are able to study the time dependence of the chemical 
rate. If we take into account the next terms in the expansion of P{t), it is easy to verify 
that we have ^^^^ < in going towards the stationary regime showing that we reach kstat 
by decreasing values of k{t) (see Appendix A). In order to to calculate the initial value k{0) 
we start from ( 13T|) . we take its time derivative and we replace d(f){t,x)/dt by its value from 
f lTU]) . After an integration by part and taking into account the boundary conditions we get 

[^^]t=o = ^</'o(O)[V.0(t,x)][,=o;t=o] = ^^o(0)[V./(a;)],=o (34) 

The last equality holds from the fact that (f){t,x) is the solution of (fTOj) verifying the initial 
condition. Since we have assumed that f{x) is entirely locahzed on the left part of x = 0, 
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we have k{0) = 0. Thus the previous results show that k{t) is a non monotonic function of 
t. 



V. CONCLUSIONS 



The main goal of this paper was to find an equation from which we may describe the 
irreversible behavior of small quantum systems. In order to do that we do not try to 
introduce a new kind of quantization as it has been proposed in the past (^) but we start 
from new and unexpected relations between dynamics and thermodynamic equilibrium. 
From the expression of the partition function we can extract a dynamics by using the 
algebraic approach of quantum mechanics or the path integral formalism. In the first case 
the dynamics is characterized by the usual evolution operator provided a rescaling of time 
is introduced while in the second approach a time irreversible partial differential equation is 
obtained. In Section 2 we have analyzed the physical content of this equation, in particular 
the relation between this equation of motion and the Schrodinger equation has been investi- 
gated. We have shown that the equation of motion describes a continuous diffusion process 
in which we keep in mind some aspects of the uncertainty relations. Although the equation 
of motion is time-irreversible it can be used to describe the thermal equilibrium and all the 
exact results concerning the thermodynamic quantities are obtained. This is based on an 
equilibrium definition that asserts that the equilibrium is established when the energy spent 
on the paths is equal to the energy put in the system from outside. From this equilibrium 
condition the thermal time r is introduced as the time on which we must explore the 
closed paths in order to have the equilibrium. This thermal time is a characteristic of the 
equilibrium. Till this point all the results investigated concern equilibrium properties. Our 
main assumption is that the equation of motion can be also used to describe the irreversible 
dynamics of systems i.e. that the equation of motion reveals the general dynamics of 
real systems. This remains an assumption since initially this equation is based on results 
associated with equilibrium states. We have established that the equation of motion is 
equivalent to a Smoluchovski equation. The density of probability, P(t,x), solution of this 
equation appears as the product of a Boltzmann distribution by a quantity counting the 
paths and depending on the initial system preparation. Nevertheless a clear meaning of 
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P{t, x) can be given. For simple systems considered here we may calculated the relaxation 
time Tr that the system needs to reach an equilibrium when it is initially in a non equilibrium 
state. We have shown that the concept of thermostat is relevant provided we have > r. 
However our results are not restricted by such a condition. We have used the Smoluchovski 
equation to investigate a standard model: the dynamics of a particle injected in a box in 
which it exists a bistable potential. We have calculated the chemical rate and we have 
shown that the chemical rate is a non-monotonic function of time for the investigated process. 



In order to conclude we may present this work on a purely deductive manner, in the 
spirit of the Feynman conjecture. Starting from a primarily discrete space-time we may 
introduce an equation of motion that retains some aspects of the uncertainty relations, this 
is an easy task since we may use the continuous limit in our level of physical description. 
This time- irreversible equation is considered as the actual equation of motion in presence 
of an external field. This equation can be used to describe exactly the thermodynamic 
equilibrium provided we consider a condition of equilibrium. This equation can be also 
used to describe irreversible situations, in this case we may transform it in a quantum 
Smoluchovski equation for which the density of probability solution of this equation has a 
clear meaning. Finally, it is possible to implement this approach by adding to the equation 
of motion an equation describing a backward motion but with the same fundametal solution. 
By forcing the system in such a way we describe a process that is reversible and our equation 
of motion is transformed into a Schrodinger equation. 



VI. APPENDIX 



In Section 4 we have restricted P{t) to Pi{t) in order to get ksta the stationary value of 
the chemical rate. To calculate the first deviation from the stationary state let consider 

P^it) = J P,,{x)cr{^^)dx + J P,,{x)c,{^^^)dx (35) 

We know (^°) that (j)o{t,x) and (j)2{t,x) are symmetric functions of x while is an- 

tisymmetric. Thus, in the domain x E [—b,a] the ratio (^;|^) is negative and (^|^) is 
positive. The constant Ci and Ci are determined by the same ratios as shown in fl2^ . We 
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can rewrite P2{t) as 

P2\t) = CLi exp 1 + 02 exp 1 (36) 

in which ai and 02 are positive numbers. The derivative of '^^dt^^ can be easily obtained 
and we get 

d -1 dP2it) 1 p^ (^1-^0) , {E2-E0). 

which is a negative quantity showing that k{t) decreases before reaching its stationary value. 
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